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The question of whether sunlight induces coherent dynamics in biological systems is under de¬ 
bate. Here we show, on the basis of an exact fully quantum mechanical treatment, that thermal light 
induces excited-state coherences in matter similar to those induced by a coherent state. We demon¬ 
strate the phenomenon on a E-type model system and a two-state Born-Oppenheimer molecular 
system. Remarkably, wavepacket-like dynamics is induced in the excited molecular potential-energy 
surface. 

PACS numbers: 33.80.-b, 42.50.-p, 42.50.Ct, 42.50.Md 


Coherence is a fundamental and central theme in quan¬ 
tum mechanics and is manifested in light-matter interac¬ 
tion through a variety of unique phenomena [T] , including 
quantum control [2] , storage and retrieval of information 
[3], and quantum state reconstruction gill]. In recent 
years, there has been a lively debate on whether sunlight- 
induced chemical reactions in biological systems evolve 
coherently and consequently exhibit extremely efficient 
energy transfer. Experimental evidence for vibrational 
and electronic coherence in such systems has been re¬ 
ported in several studies [SHS]. However, it has been ar¬ 
gued that since such studies employ coherent laser pulses 
rather than sunlight (so-called thermal or chaotic light) 
to initiate material dynamics, the results obtained can¬ 
not establish conclusively that sunlight-induced biologi¬ 
cal processes evolve coherently [14]. 

Theoretical studies of the excitation of atomic and 
molecular systems by thermal light have reached contra¬ 
dictory conclusions as to whether the induced dynamics 
is coherent. Most of these studies are approximate in 
some sense, as they involve either semiclassical or per¬ 
turbative treatments, and some of them rely on the pres¬ 
ence of an intrinsic coupling between the material excited 
states, such as vacuum- or decay-induced coherence [15]- 
[n\ . Others, treat the material system only following an 
assumed coherent excitation [53]. None of these studies 
employ an exact and fully quantum mechanical treatment 
of the interaction of thermal light with matter without 
such limiting constraints. Interestingly, the entanglement 
(correlations) between separated material sites has also 
been studied in this context [M] [^^ . 

Since the interaction of thermal light with matter is 
a fundamental issue of enormous potential implications, 
the need for a definite conclusion on whether thermal 
light can induce coherent material dynamics is of pri¬ 
mary importance and significance. In order to reach a 
definite concfusion, we treat the problem starting from 
basic considerations which, to the best of our knowledge, 
has never been done before. 

We employ an exact (non-perturbative) dynamical 
treatment for the interaction of a material system 


(atomic and molecular) with a thermal state of the elec¬ 
tromagnetic field, where both the material and light sub¬ 
systems are quantized. Concretely, we treat the light- 
matter system as a bipartite composite system whose 
subsystems interact via a Jaynes-Cummings-type |26j in¬ 
teraction model. We thus express the initial system by 
its density matrix, p(0) = Pf{0) C Pm(0) (where F and 
M stand for ‘field’ and ‘material’, respectively), and fully 
propagate it in time according to p{t) = U(t)p{0)W (t), 
where U{t) is the system propagator. In the case we con¬ 
sider, the material subsystem is initially in its ground 
state, while the field is in a single-mode statistical mix¬ 
ture of Eock states. 


As we describe below in detail, we apply the interaction 
model to a three-level E-type system as a basic model, 
and, for the first time to the best of our knowledge, 
to a two-state Born-Oppenheimer molecular system for 
which we propagate the full density matrix, represented 
in the electronic C bond-coordinate C photon-Eock prod¬ 
uct space. To conform with the conditions of a natural 
process in our simulations, we consider the weak-field 
regime and initial average photon number characteristic 
of sunlight (about 500nm wavelength at 6000 K, [57]). 
We find that thermal light induces excited-state coher¬ 
ent dynamics in both the atomic and molecular subsys¬ 
tems, as it also does in the regime of strong field and 
high average photon number. We compare this result 
with that obtained for the interaction of the material 
systems with a coherent-state. Remarkably, we find that 
both scenarios exhibit almost identical excited-state co¬ 
herence features. On the other hand, in contrast to a 
coherent-state, thermal light does not induce coherence 
between the material ground and excited states. 


As a prototype for a molecular system, we consider a 
three-level E-type model system interacting with a single 
mode of the radiation field. The Jaynes-Cummings-type 
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system Hamiltonian is 

H = ^a;i|*)(j|+w(a’^a+i) 

i 

+ ^ Ai (|«)(5|a+l5)(*|a^) > (1) 

i¥^a 

where we have invoked the dipole and rotating-wave ap¬ 
proximations, and set h = 1. The index i stands for 
the ground and two excited states of the V system, g, e 
and /, respectively. The operators and a are the pho¬ 
ton creation and annihilation operators, respectively, op¬ 
erating on the photon states {!«-)}. The parameter Ai 
is the matter-field interaction constant. Note that this 
model excludes any intrinsic coupling between the ex¬ 
cited states. 

The material system is initially in its ground state 
/Om(0) = \g){g\- For the field we consider two possible 
initial states: a coherent state and a thermal state. A co¬ 
herent state takes the form |a) = e“ 2 l“l 
where jap = h is the field average photon number. The 
corresponding density matrix is p_f( 0) = |a)(a|. For the 
thermal state pf(0) = J2nPn\n'){M^ with p„ = ■p:;|^s+T, 
and the average photon number h, given by the Boltz¬ 
mann distribution m- Thus, the initial state of the 
composite system is p(0) = pf(0) 0 Pm{0), and the state 
of the system at any time in its evolution is given by 

p{t + At) = (2) 

The state of each subsystem is obtained by tracing p{t) 
over the space of the other subsystem; that is, 

PA{t) = Trs [p(t)]. (3) 

In the supplemental material, we give an analytical ex¬ 
pression for p{t) for the V system, as well as for the co¬ 
herence element of the material subsystem for the fully 
resonant case where ui — = 0. Furthermore, we com¬ 

pare the analytical and numerical calculations in order 
to justify the numerical simulations presented here. 

We now compare the V system dynamics induced by 
a coherent state with that induced by a thermal state. 
For both scenarios, the energy separation between the 
|e) and |/) states is set to 250 cm“^, and the field fre¬ 
quency is tuned between them. We set \ej = 10“^, and 
h = 0.008 (obtained for 495.9 nm and temperature 6000 
K for the field mode). We construct the initial states 
of the light with the first ten Fock states. In Fig. [l] 
we show the material excited states populations (top 
panel), the real and imaginary parts of the material ex¬ 
cited states coherence element, pMje(t) (middle panel), 
and the trace of the material (and field) reduced den¬ 
sity matrix, Tr[pM,F(i)]> of if® square, Tr \_p\[ p it)] 
(bottom panel), as obtained for the interaction with the 
coherent state. In Fig. we present the same dynam¬ 
ical measures for the interaction of the V system with 


thermal light; in the lower panel, however, we show only 
Tr [pIi (t)] and Tr [p|’(t)] which are generally different 
for a mixed state. 




t(fs) 

FIG. 1: (color online). Populations (top), coherences (middle) 
and partial traces (bottom) obtained for the interaction of a coher¬ 
ent state with a three-level F-type system. 


Figures [T] and show striking similarity between the 
material populations and, most remarkably, the excited 
state coherence for the coherent and thermal light sce¬ 
narios. In this context, it is interesting to refer to the 
similarity in the interaction of a coherent and thermal 
states with a two-level system as noted by Cummings 
[25] . In contrast, Tr Tr [p|n(t)] show dif¬ 

ferent behavior. In the mixed state scenario (interaction 
with the thermal light), the two oscillate out-of-phase (as 
if purity is exchanged between the subsystems), while in 
the pure scenario (interaction with the coherent state) 
both are identical, as is well known for a pure system. 
We note that throughout the interaction of the thermal 
light with the V system, pF{t) remains diagonal as we 
show analytically in the supplemental material. 

Another major difference between thermal and coher¬ 
ent state of the light is that while the latter induces co¬ 
herence between the material ground and excited states, 
the former does not. (We derive this result analytically in 
the supplemental material). We explain this phenomenon 
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FIG. 2: (color online). Populations (top), coherences (middle) 
and partial traces (bottom) obtained for the interaction of thermal 
light with a three-level F-type system. 


in the following way. According to our model, the inter¬ 
action of a single Fock state, say |n), with the mate¬ 
rial ground-state \g) produces the general superposition 
a\g,n) + (5|e) -I- c|/))|n — 1), which in the material sub¬ 
space (that is, after tracing over the field subspace) cor¬ 
responds to a mixed state of the material ground and ex¬ 
cited states. However, the material excited-states are in 
a coherent superposition with coherences between them. 
Thermal light is a statistical mixture of single Fock states 
and, therefore, induces only excited-state coherences in 
the material subsystem. The absence of coherence be¬ 
tween the ground and excited states is related to the 
vanishing of the material dipole moment in this scenario 
as pointed out in [35]. In contrast, since a coherent state 
is a coherent superposition of Fock states, its interaction 
with the material ground-state will generate coherence 
between the material ground and excited states (since 
they share common photon states), and also among the 
excited states themselves. Note that a Schrodinger cat 
(coherent) state [37j is expected to behave just like a 
thermal state, as noted elsewhere [30j . 

The material excited-state coherence, induced by the 
thermal light, is obtained through an exact quantum 


mechanical treatment without any intrinsic coupling be¬ 
tween the excited states. Such couplings were introduced 
in a number of previous studies, and were necessary to 
allow the material excited-state coherence. Our exact 
treatment proves such coherence is obtained naturally 
via the interaction with thermal light, similarly to that 
obtained in the interaction with a coherent state. In the 
supplemental material we show that similar dynamical 
characteristics are obtained for the interaction of a the 
V system with two-mode thermal (and coherent) state. 

A realistic molecular system is by far more significant 
for understanding matter interaction with thermal light. 
We therefore apply a similar fully exact quantum me¬ 
chanical treatment to a molecular system of two Born- 
Oppenheimer potential energy surfaces (PESs). Specif¬ 
ically, the field is represented in Fock space, and the 
molecular system is spanned by the two-dimensional elec¬ 
tronic space and the continuous coordinate space for the 
nuclei separation. Thus the state of the system is given 
by the tensor product of these three spaces. 

The Hamiltonian of the entire system is 

Hmoi= Hi{r)\i){i\ + uj{a'a + |) 

i^g,e 

+X{\e){g\a+\g){e\a^) , (4) 

where Hi{r) = Vi{r) + T is the nuclear Hamiltonian of 
the ground and excited electronic states (I 5 ) and |e)), 
with the PES Vi{r) and the kinetic energy operator 
T = parameter m is the molecular reduced 

mass. We denote the (vibrational) eigenstates of Hi{r) 
by where v is the vibrational quantum number. 

In writing the Hamiltonian of Eq. |^we invoke the dipole, 
Condon, and rotating-wave approximations. 

We examine the molecular dynamics induced by 
excitation with both thermal and coherent states, 
as before. The initial molecular (ground) state is 
/Om( 0 ) = Ig,' 0 g,o )(57 ^g,o|; and the initial state of the 
field pf( 0) is either a coherent or thermal state, as 
described above. The initial state of the system is 
p(0) = pf(0) 0 Pm(0), and its state at any time in the 
course of interaction is 

pit + At) = (5) 

It is convenient to employ the “split operator” method 
m in propagating p, and split the kinetic en¬ 
ergy term of the Hamiltonian from the rest of the 
terms. Thus, for each time-step the propagator is 
^ where Hmoi is the full 

Hamiltonian of Eq. without T. In practice, we sim¬ 
ulate the interaction of thermal light with a two-state 
Born-Oppenheimer model system of two one-dimensional 
Morse-type PESs, Hi(r) = Di{\ — 

potentials parameters, based on these of the Li 2 molecule, 
are given in the supplemental material). The frequency 
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of the exciting field is tuned between the third and fourth 
vibrational eigenstates of the excited state. The initial 
state pf(0) = Yl,nPnW){''^\ is constructed with the first 
seven Fock states. We set A = 10“"^, and n = 0.0072 
(obtained for 486.1 nm at 6000 K for the field mode). 

In Fig. we show snapshots of the excited-state 
material reduced density matrix, pM,ee{t)- In the left 
column we show pM,ee{t) in the coordinate representa¬ 
tion, while in the right column we show its projection 
onto the first 15 excited vibrational eigenstates. The 


in-5 



FIG. 3: (color online). Snapshots of PM,ee(i) for thermal light 
interacting with two-state Born-Oppenheimer molecule. The left 
column shows PM,ee(f) in coordinate space (with the time specified 
in femtoseconds), and the right column shows its projection onto 
the vibrational eigenstates of the molecular excited-state. 


clear observation that emerges from these snapshots is 
that thermal light induces excited-state vibrational co¬ 
herent dynamics. An almost identical coherence fea¬ 


ture is obtained by the interaction of a coherent state 
with the molecular system (not shown here). The spa¬ 
tial representation indicates wavepacket-like dynamics, 
unlike that reported recently on the basis of a semi- 
classical treatment [5D]. (When the field is resonant 
with a vibrational eigenstate of the excited-state, the 
state created, both by the interaction with coherent and 
thermal states, is the corresponding vibrational eigen¬ 
states). The difference between the interaction with 
coherent and thermal state is that the former also in¬ 
duces coherence between the ground and excited elec¬ 
tronic states, while the latter does not, as discussed above 
for the V system. To have a quantitative measure of the 
wavepacket-like dynamics, we calculate the excited-state 
correlation function C{t) = Tr [pM,gg{0)pM,eeit)] and its 
Fourier-transform < 7 ( 0 ;), shown in the supplemental ma¬ 
terial. The spectral signature is characterized by the har¬ 
monics of the energy difference between the excited-state 
vibrational eigenstates ojex ~ 0.0011 a.u. 

The excited-state molecular coherent dynamics, ob¬ 
tained in our treatment of the interaction with thermal 
light contradicts previous studies that either excluded co¬ 
herence in the material system, or allowed the existence 
of coupling between the excited states that was neces¬ 
sary for the generation of coherence mm]- According 
to our treatment, the generation of coherence is an in¬ 
trinsic property of the interaction between the thermal 
light and material system. 

In conclusion, the question of whether sunlight induces 
coherent dynamics in biological systems has been contro¬ 
versial for many years now, with contradicting answers. 
Semiclassical perturbative studies excluded the possibil¬ 
ity of thermal light-induced coherence, while other stud¬ 
ies found it possible, provided there is an intrinsic source 
that couples the excited states of the system. 

Here we have shown, on the basis of an exact fully 
quantum mechanical treatment with a Jaynes-Cummings 
type of interaction model, and under conditions that con¬ 
form with sunlight-induced natural processes, that ther¬ 
mal light does indeed induce excited state coherent dy¬ 
namics, both in a three-level H-type system and in a 
two-state Born-Oppenheimer molecular systems. The 
induced coherence is similar to that induced by a co¬ 
herent state. Remarkably, the thermal light is shown to 
induce wavepacket-like dynamics in the molecular sys¬ 
tem. Unlike a coherent state, however, the thermal light 
does not induce coherence between the material ground 
and excited states. Most significantly, and in contrast 
to previous studies, the generation of the excited-state 
coherences in our model requires no intrinsic coupling 
between the excited states. In general, thermal light ex¬ 
cludes coherence between the states it couples directly, 
and permits coherence between material states which are 
not coupled directly by the field or by any other coupling 
mechanism. Interestingly, a Schrddinger cat (coherent) 
state is expected to induce coherences similar to thermal 
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light. On the basis of recent studies EH, we expect that 
the interaction of the system with environment would not 
eliminate the induced coherence. In addition, based on 
our two-mode calculations, we expect that a continuous¬ 
mode of thermal light, interacting with matter, should 
induce similar type of coherences. 

The results we present in this Letter may suggest a 
signihcant role that a statistical source of modes can 
have in coherent dynamics, as was pointed out in a 
similar context [23]. Consequently, our results may 
have direct impact on the understanding the dynamical 
characteristics of sunlight-induced biological processes. 

^Author to whom correspondence should be addressed, 
david. avisar@gmail. com 
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SUPPLEMENTARY MATERIAL 

Analytical dynamical expressions for the interaction 
of the V system 

In the interaction picture, the Hamiltonian for the V- 
type three-level system, assuming uj — ujej = 0 and a 
common parameter A, is given by 

Hi = (l*)(5|a+ l5)(i|a^) 

i=ej 

( 6 ) 

Employing Taylor expansion, we get the following ma¬ 
trix for the propagator in the three-dimensional 

material space 


(Cf{t) Cfeit) Sfgit)\ 

= 4/(f) Ceit) Seg{t) , 

\Sgfit) Sg^it) Cg{t)j 


with the matrix elements given by 


Cfit) 

Cfe{t) 

Sfg{t) 

Sgfit) 

Cgit) 


= Ce{t) = ^ cos(V2aalAt) -I-1 
= Cef{t) = ^ cos{V2aa'^ Xt) — 1 


= Seg{t) = 


= Sgeit) = 


—i sin (•\/2aal At), 
\/2 V afit 
—i sin(-\/2alaAt) 


V2 

= cos('\/2alaAt). 


(7) 


( 8 ) 


where we have omitted, for convenience, the explicit time 
dependence of the matrix elements of Eq. 

First we derive an expression for the coherence element 
between the ground and excited state for the case of a co¬ 
herent state of the field. We then show that this element 
is identically zero for the thermal state of the field. Let 
us consider the element —CgPF{0)Sge- For the coherent 
state of the field we plug pf( 0) = |Q!)(a| into this element 
and trace over the field space. As a result, we get 


PM,gf{t) 


Ttf \^-Cg(t)\a){a\Sgeit) 
-^{m\Cg(t)\a){a\Sgeit)\m) 


-{a\Sgeit)Cg{t)\a) 



E 


*)«+! ^ 
— , =sin 

y^n\{n + 1)! 


\/2(n -I- l)At 


xcos(v^At). 


( 11 ) 


If we plug into —CgPFiO)Sge the expression for the ther¬ 
mal state of the field, trace over the 

field space, the coherence element pM,gf{t) is identically 
zero (as well as for the eg element): 


PM,gf{t) 


Ti-f 


-Cgit)'^Pu\n){n\Sge{t) 

n 


-'^Pn{m\Cg{t)\n){n\Sge{t)\m) 

n,m 


-'^Pn{n\Sge{t)Cg(t)\n) = 0 (12) 

n 


The coherence element pMjeit) is obtained by tracing 
over the field subspace of the matrix element pijf,(t): 


The operations of each of the propagator matrix elements 
on Fock states are then 

Cfit)\n) = i cos(\/2(n -|- l)At) -I- 1 |n) = C'/,„(t)|n) 
Cfe{t)\n) = i cos(\/2(n -|- l)At) - 1 \n) = Cfe,n{t)\n) 
Sfgit)\n) = ^sin(v^At)|n - 1) = S'/g,„(t)|n - 1) 

Sgfit)\n) = ^sin(\/2(n -h l)At)|n -f 1) = S'g/,„(t)|n -h 1) 
= cos{V^Xt)\n) = Cg_„(t)|n). 

(9) 


PMJe{t) = TrF ^-Sfg{t)pFiO)Sgeit) 

= -J2{rn\Sfg{t)pF{0)Sge{t)\m). (13) 

m 

For a coherent state, we plug Pf( 0) = |Q:)(a| (with 
la) = e“ 2 l“l 1^9 1^ set the fol¬ 

lowing expression for the excited-state coherence element 


PM,fe{t) = -e ^-^^sin^(v^At) 


= -e —-sm^{V^Xt). (14) 

2 nl 


For the initial state p(0) = \g){g\ ® Pf( 0) of the com¬ 
posite system, the density matrix at any time t of the 
evolution, in the interaction picture, is 

(~^fgPF{^)Sge —SfgPF{0)Sge SfgPFiO)Cg\ 

= I —SfgPF{0)Sge —SfgPF{0)Sge SfgPF{0)Cg I , 

\—CgPF{0)Sge —CgPF{0)Sge CgPF{0)Cg/ 


For the case of thermal field, with pf( 0) = X)„P„|n)(n|, 
we get 


Pm 


Je{t) = - y^PnSin^(v^At) 


( 10 ) 


- E 


(1 -h h)”+i 


sin^ (V^Xt). (15) 
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In addition, we can also show that during the interac¬ 
tion of the thermal light with the V system, the reduced 
density matrix of the light, ppit), remains diagonal. To 
show that we trace pi{t) of Eq. 10 over the material sub¬ 
space, and use the the equalities of Eq. The result 
obtained is 


PF{t) = TiM[pi{t)]=Y^{i\pi{t)\i) 

i 

= -‘^Sfg{t)pF{0)SgS) + Cg{t)pF{0)Cg{t) 


= ^p„sin^(v^Af)|n — l)(n — 1| 

n 

+ y^PnCos^(v^Af)|n}(n|, (16) 

n 


which represents a diagonal reduced density matrix 

PF{t). 

We use the above analytical expressions for the ex¬ 
cited states coherence to compare with our numerical re¬ 
sults and verify our numerical simulations for the den¬ 
sity matrix dynamics. In Eig. we show the real part 
of the excited state coherence element of the density ma¬ 
trix PMjeit), obtained for the field initially in a coherent 
state. In blue is the result obtained from the numerical 
propagation, whereas in red is the result obtained using 
Eq. 14 In Eig.([^, we show the comparison of the an- 



FIG. 5; (color online.) Real part of pM,fe{t) obtained numer¬ 
ically, after Eq. 2 of the main text (blue) and analytically, after 
Eq. |15| (red) for the interaction of the thermal light with a V sys¬ 
tem. 


Interaction of V system with two-mode held. 

Here we show that, similar to the interaction of a 
singl-mode states of the light, the interaction with two¬ 
mode states of a coherent state and thermal state, induce 
excited-state coherence. 
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FIG. 4: (color online.) Real part of pMJei^) obtained numer¬ 
ically, after Eq. 2 of the main text (blue) and analytically, after 
Eq. |14| (red) for the interaction of a coherent state with a V sys¬ 
tem. 
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alytical (following Eq. 151 and numerical calculation for 
the excited state coherence obtained for the thermal field. 
Note the remarkable similarity of the results obtained for 
the coherent state of the light and the thermal state! 
Figures @ and (§ confirm, unequivocally our numerical 
propagation for the density matrix. 


FIG. 6: (color online.) Populations (top), coherences (middle) 
and partial traces (bottom) obtained for the interaction of a two¬ 
mode coherent state with V system. 


In Fig.j^we show the dynamical measures, as discussed 













































Spectrum of the excited-state coherent dynamics 


in the main text, of the interaction of a system with 
two-mode of a coherent state. In Fig. we show the same 
dynamical measures obtained for the interaction of a V 
system with two-mode thermal state of the light 



FIG. 7: (color online.) Populations (top), coherences (middle) 
and partial traces (bottom) obtained for the interaction of a two¬ 
mode theraml light with V system. 


In Fig. we show the spectrum obtained by Fourier- 
transform of the excited-state correlation function 
C{t) = Ti [pM,gg{0) Pm, ee{t)]- The spectrum reveals the 
harmonics of the energy difference between the excited- 
state vibrational eigenstates ujex ~ 0.0011 a.u. 



FIG. 8: (color online). The spectrum obtained by Fourier- 

transform of the excited-state molecular auto-correlation function 
C(i) = Tr [pM,gg{0)pM,se{t)]- 


Potentials parameters 


Below are the potentials parameters used for the two- 
state Born-Oppenheimer molecular dynamics. The PESs 
are essentially the X and A states of the Li 2 molecule, 
but only with the A state shifted in energy |T]. 


TABLE I: The parameters, in atomic units, for PESs used for the 
molecular simulations. 



Vg 

K 

D 

0.0378492 

0.0426108 

b 

0.4730844 

0.3175063 

ro 

5.0493478 

5.8713786 

T 

0 

0.0911267 


[1] G. Herzberg, in Molecular Spectra and Molecular Struc¬ 
ture; I. Spectra of Diatomic Molecules, Krieger Publishing 
Gompany, Malabar, Florida, (1950). 









































